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Abstract 

The precision measurement of the hyperfine splitting Aiip{lP, cc) = McogiXcj) — M{hc) = 
—0.5 lb 0.4 MeV in the Fermilab-E835 experiment allows to determine the gluonic condensate G2 
with high accuracy if the gluonic correlation length Tg is fixed. In our calculations the negative value 
of Ahf = -0.3 ± 0.4 MeV is obtained only if the relatively smah Tg = 0.16 fm and G2 = 0.065(3) 
GeV^ are taken. These values correspond to the "physical" string tension (a ~ 0.18 GeV^). For 
Tg > 0.2 fm the hyperfine splitting is positive and grows for increasing Tg. In particular for Tg = 0.2 
fm and G2 = 0.041(2) GeV^ the splitting Ahf = 1-4(2) MeV is obtained, which is in accord with 
the recent CLEO result. 

PACS numbers: ll.15.Tk, 12.38.Lg, 14.40.Gx 
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I. INTRODUCTION 



Recently, the CLEO Collaboration has presented preliminary results on the mass of the 
hc{l^Pi) resonance in the reaction tp' vr^/ic Q]. It appeared that the central value of 
the observed mass, M{hc) = 3524.4 ± 0.9 MeV, is about 1.5 MeV lower than the one 
found in the recent Fermilab-E835 experiment, in the pp he rjc'y 87 reactions, 
where M{hc) = 3525.8 ± 0.4 MeV is measured jjj. Consequently, the central values of the 
hyperfine (HF) splittings, Ahf(1-P, cc) = McogiXcj) — M{hc), corresponding to the measured 
masses M{hc), have different signs: 

Aij^^(exp) = -0.5 ±0.4 MeV (E835), (1.1) 
AgJ(exp) = +0.96±1.03MeV (CLEO), (1.2) 

although they are compatible within 2a. (In Eqs. fll.lU1.2|) the value Mcog(Xcj) = 3525.32 ± 
0.13 MeV is taken from PDG Js]). The small numerical difference between these experimental 
values makes an essential difference for the theoretical interpretation, as we will show in this 
paper. 

Already many years ago it was understood that the sign and small magnitude of Ahf(I-P) 
in charmonium occur due to the almost equal (and small) contributions from the perturbative 
(P) and nonperturbative (NP) HF interactions 

Ahf(1P) = A^p(1P) + A^P(1P), (1.3) 

where the perturbative contribution is always negative 0, 0]. Just due to the cancella- 
tion between the negative perturbative contribution and the positive NP one, the value 
Ahf(1P, cc) = —1 MeV has been calculated in [6] for a value of the gluonic condensate 
G2 = 0.042 GeV*^, while in Ref. the same HF splitting has been obtained for the signif- 
icantly smaller value G2 ~ 0.02 GeV^. The reason behind this difference will be explained 
below and comes from the fact that the gluonic condensate actually enters AyF(lP) in the 
combination ^2^3 {Tg is the gluonic correlation length) and therefore the extracted value 
of G2 depends also on the correlation length Tg used. Unfortunately, at present there is a 
large uncertainty in the value of the gluonic condensate, even in the framework of the same 
approach, like QCD sum rules. Values ranging from G2 = 0.012 GeV'* up to G2 = 0.07(3) 



GeV*^ are used 
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The new experimental value of the he mass j2| (measured with an accuracy of 0.4 MeV) 
allows one to determine G2 with better precision. In this paper we show why a precise knowl- 
edge of Ahf(I-P) in charmonium is so important for a fundamental theory and explicitly 
extract the gluonic condensate G2 from existing experimental data. 



II. THE PERTURBATIVE HF INTERACTION 



We consider here the perturbative HF interaction in one-loop approximation which is 
well known. The splitting Ahf(I-P) is given by the expression jsl: 



4 3 



{r'')iP. (2.1) 



It is important to take the HF splitting just in the form 1)2.11) while another (approximate) 
definition suggested in Ref. P], 



1 _ !^ 
4 3 



TT 



{M(3p2)-M(3Po)}, (2.2) 



cannot be used because it is valid only in lowest order and neglects second order corrections 
(~ «s/7i") and the NP contribution. It will be shown later that the contribution of the 
neglected terms to the difference M{xc2) — M{xco) = 3a + 0.9c (where a and c are the 
spin-orbit and the tensor splitting) is about 26% and has negative sign (see the numbers in 
Table P). 

It is also important that if relativistic corrections are taken into account, the current (pole) 



quark mass in 1)2.1)) must be replaced by the average kinetic energy, ujg = (yP^ + ^g)ip, 
(usually called the constituent quark mass). This modification can be rigorously derived 
within the Fock-Feynman-Schwinger representation of the gauge-invariant meson Green 
function, when the spin-dependent interaction can be considered as a perturbation jic| . 
Then instead of ()2.1)) one must use 

ASF(lP) = -|^^(r-^)iP (%=4). (2.3) 

As seen from Eq. fl2.3)) the HF splitting strongly depends on the coupling as(yu), where the 
scale /i in Q;s(/i) cannot be arbitrary. It is clear that for the ^Pi state the scale should be 
the same as in the coupling apsifJ') used in the fine structure (FS) splittings of the Xc(^-Pj) 
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mesons. Fortunately, apsif^) can be directly extracted from the experimental values of the 
spin-orbit and tensor splittings: as derived in the following relation is valid 

2 / N vr^g{r?c(exp) - |aNp(lP)|} 
-.M = ^J^^ • (2.4) 

Here ?7c(exp) = |c(exp) -a(exp) = 0.024(1) GeV; the tensor splitting c(exp) = 0.039(1) GeV 
and the spin-orbit splitting a(exp) = 0.0346(2) GeV are calculated from the Xcj masses. 
In Eq. fj2.4p aNp(l-P) is the Thomas precession (NP) term: 

aNp(lP) = --^(r-i)iP, (2.5) 

while /4 is the matrix element (m.e.) entering the second order perturbative part of the 
parameter rjp: 



n - ^r n - ^ ^(2) (2) _ 2a|s(^ 

rjp — -cp — ap — -Cp — ttp — — /4. [Z.io) 

2 2 Tiu^ 

For the extraction of the strong coupling constant it is essential that the m.e. /4 does not 

explicitly depend on the renormalization scale /i: 

/4 = 1.97834(r-3)ip - {r'^ ln(mr))ip {rij = 4). (2.7) 

The derivation of the HF interaction in Ref. shows that the mass m{m), entering /4, is 
the current (Lagrangian) mass (m = 1.16 GeV and rric (pole)=1.45 GeV are taken here). 

The analysis in Ref. has shown that the extracted value of q;fs(/^fs) (in the MS 
scheme) and the scale /ips significantly differ if the constituent mass Uc (about 1.6 ^ 1.7 
GeV) instead of the pole mass (about 1.4 -^1.48 GeV) enters the expression (|2.4j) . For 
example, the precise description of the FS splitting (with an accuracy better than 1%) for 
mc (pole) = 1.45 GeV in Eq. gives 

apsif^i) = 0.358 with /xi = 0.70 GeV. (2.8) 



In relativistic calculations with ci;c(lP) = (yp^ + ^c)ip — 1-66 GeV the value of the cou- 
pling constant apsifJ'R) extracted from Eq. ()2.4|) appears to be significantly larger, 

«Fs(/^ii) = 0.514; fiR = 0.51 GeV. (2.9) 

From our point of view it is the second choice which should be prefered, because the value 
ct^Fs(/^i) in Eq. ()2.8p is too small for such a small scale as /xi = 0.70 GeV. It is known that 



TABLE I: The spin-orbit and tensor splittings (in MeV) of the Xcj mesons (the static potential 



131] with the parameters mc(pole)= 1.45 GeV, ujc = 1.66 GeV, a = 0.18 GeV^, 



is taken from 
A^(2-loop)= 267 MeV, aFs(/^FS = 0.51 GeV)=0.514. 





1st order term 


2d order term 


NP term 


a(tot)= 34.58 
a(exp) = 34.56(19) 


Op^ = 51.92 


af'^ = -4.24 


«NP = —13.10 


c(tot)= 39.12 
c(exp) = 39.12(62) 


c}^ = 34.61 


= 4.51 


CNP < 1 MeV") 



°^ The reasons why cnp(I-P) is smah, are discussed in 
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3- 



the strong coupling in the MS scheme is already rather large at the scale = 1.777 GeV 
(asiMr) > 0.33 pI), while at the scale /ips = 1-0 GeV the coupling aFs(1.0 GeV)^ 0.40 
has been obtained in jl^ in the analysis of FS splittings of the 2P state in bottomonium. 
Therefore the values given in Eq. ()2.9|) will be taken here. They result in the following FS 
splittings of the Xcj mesons (with an accuracy better than 1%): 



c(lP) = 39.12MeV, Cexp(lP) = 39.12 ± 0.62MeV, 
a(lP) = 34.58MeV, acxp(lP) = 34.56 ± 0.19MeV. 



(2.10) 



The first- and second-order perturbative and NP terms in the FS splittings: a (IP) = a^^-* + 



,(2) 



+ onp, c(1P) = c!p^ + c!p^ (cnp is very small 1^), are given in Table HI 



Having obtained a precise description of the FS of the Xcj mesons, we can expect that 
the HF splitting (for the same set of physical parameters) is also determined with good 
accuracy. From the expression ()2.3|1 (with Uc = 1.66 GeV, anp^fiR) = aFs(/^i? = 0.51 
GeV)= 0.514, {r~^) = 0.139) it follows that in charmonium the perturbative contribution 
has the value 

APp(lP) = -4.1 MeV, (2.11) 

which is five times larger than the one obtained in the experiments Cm the CLEO 

experiment |3| Ahf is positive, 1.0 ±1.0 MeV). 

The matrix elements Eq. ()2.7|1 are calculated here utilizing the solutions of the spinless 
Salpeter equation with a static potential-linear plus gluon-exchange term, where in two-loop 
vector coupling the asymptotic freedom behavior at small distances and the freezing of the 
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coupling at large distances are taken into account. The most important matrix elements are 
(r-i) = 0.405 GeV, {r-^)ip = 0.139 GeV^, and {r-'^ \n{mr)) = 0.095 GeV^. Note that the 
m.e.(r~^) in the relativistic case is about 30% larger than in nonrelativistic calculations. 



III. THE NONPERTURBATIVE HF INTERACTION 

Spin-dependent NP potentials have been introduced in [ij, Q. With the use of the 
vacuum correlation function (v.c.f.) D{x) the HF interaction is written as 

1 



with 



NP 
HF 



6 



vrir) 



3uj 



2 M 



duDlV + z/^ 



(3.1) 



(3.2) 



The contribution of the other correlator, Di{x), has been neglected in Eq. ()3.2p . because 
in the unquenched case Di{x) is small, even compatible with zero jl3]. The v.c.f. D[x) 



is shown to behave as an exponential: D{x) = dexp{—x/Tg) at a; > 0.2 fm |l2| {Tg is 



the gluonic correlation length), while at smaller x < tq {tq < 0.2 fm) D{x) should have a 



SD{x) 



< and ^ 



x=0 



0, established in 



)lateau to satisfy the necessary conditions 

The value of D{x) at the origin, -D(O) is related to the gluonic condensate G2 



IT 



(3.3) 



and we assume that D{ro) = -D(O), so the factor d in front of the exponent is d 
D{0) exp{ro/Tg). At the origin the HF interaction Eq. ()3.1|) has the value ^[^(r = 0) 
n^G2 (ro + Tg)/(9a;2). 

The string tension in the confining potential is defined through the same v.c.f. D{x): 

poo roo 

a = 2 dX duD{Vx 
Jo Jo 

and for the form adopted for D{x) we find 



^2^Z/2) 



(3.4) 



i + ^ + if^ 



2 \Tn 



(3.5) 



^9 ^ x-^a/ 

From this relation and taking the "physical" value a ~ 0.178(8) GeV^ one can determine 
the gluonic condensate for different values of Tg. However, the gluonic correlation length is 



not known with good accuracy and at present different values (between 0.16 fm and 0.26 
fm) have been obtained in lattice QCD 17 1 and in the vacuum correlator method jioj ]. 
For comparison we use here two values, Tg = 0.16 fm from Ref. [w| and 0.2 fm. We do not 
consider here the value Tg ^ 0.3 fm, because our analysis shows that for such a value of Tg 
the NP contribution to the HF splitting is too large (A2p(lP) ^ 10 MeV) so that the total 
HF splitting is positive, about 6 MeV, i.e. several times larger than the experimental value 
given in Eq. p.2|) . 

Then from the relation dSUl) with a = 0.178(8) GeV^ it follows that 

G2 = 0.065(3) GeV^ (ro = Tg = 0.16 fm), 

G2 = 0.041(2) GeV^ (ro = Tg = 0.20 fm). (3.6) 

The HF splitting, calculated for the interaction ()3.2|) . reduces to the expression 0, 

Ahf(I^) = y §(^0 + Tg)J, J = (rK, (^^) (3.7) 

The accuracy of this approximation is better than 5%. The matrix element J, J{Tg = 
0.16 fm) = 0.092 GeV"^ and J{Tg = 0.20 fm) = 0.17 GeV-\ strongly depends on Tg. Then 
taking G2 from (j3.6|l one obtains 

A^J^(IP) = 3.8(4) MeV {Tg = 0.16 fm), 

A^JP(IP) = 5.5(3) MeV {Tg = 0.20 fm). (3.8) 

Thus the magnitude of the NP contribution ()3.8|) appears to be larger (smaller) than the 
perturbative term ()2.11|) for larger (smaller) gluonic correlation length. Therefore the total 
HF splitting Eq. p.3|l has different signs for Tg = 0.16 fm and Tg = 0.2 fm, 

Ahf(IP) = -0.3(4) MeV {Tg = 0.16 fm) (3.9) 
Ahf(IP) = +1.4(2) MeV {Tg = 0.20 fm). (3.10) 

It is amusing to notice that the HF splitting for Tg = 0.16 fm exactly coincides with the 
experimental value obtained in the E835 experiment |2j, while the positive splitting ()3.10p 
is close to the value obtained in the CLEO experiment Q|. To distinguish between these 
two possibilities Ahf(IP) needs to be measured in charmonium with a better accuracy. 
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IV. CONCLUSIONS 



Our analysis has shown that the "physical" value of the string tension cannot unam- 
bigously fix the gluonic condensate and only the product G2Tg can be determined. 

The HF splitting between the c.o.g of the Xcj niesons and /ic(^-Pi) provides an additional 
opportunity to extract the gluonic condensate with good accuracy, due to the cancellation 
between the negative perturbative and positive NP contributions which both have small 
magnitude. 

To calculate the HF splitting for the IP states in charmonium we use exactly the same 
coupling as{fi) and matrix elements as in the fine structure analysis of the Xcj niesons where 
a high accuracy (< 1%) is reached. Therefore we estimate the accuracy of our calculations 
of Ahf(I-P) to be equal to 0.3 MeV (0.4 MeV) for the gluonic correlation length Tg = 0.16 
fm {Tg = 0.2 fm). 

An additional restriction is also put on the gluonic condesate-it should correspond to the 
"physical" string tension, cr = 0.18 GeV^, used in the static potential. 

Then a negative central value of Ahf(I-P) = —0.3 MeV, as in the E835 experiment 
P], is obtained for G2 = 0.065 GeV^(Tg = 0.16 fm). For the larger correlation length, 
Tg = 0.2 fm, and the smaller value G2 = 0.041 GeV^ the HF splitting appears to be 
positive, Ahf(I-P) = +1-4 MeV. So it is of prime importance for the deduction of the 
values of the gluonic condensate and correlation length that the discrepancy between the 
two experimental results be removed. 
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